arXiv: 1506.00205v2 [cond-mat.stat-mech] 21 Nov 2015 


Spontaneous Symmetry Breaking and Phase Coexistence in Two-Color Networks 


V. Avetisov 1,2 , A. Gorsky 3,4 , S. Nechaev 5,6 , O. Valba 1,2 
1 N.N. Semenov Institute of Chemical Physics of the Russian Academy of Sciences, 119991, Moscow, Russia, 

2 Department of Applied Mathematics, National Research University Higher School of Economics, 101000, Moscow, Russia, 
3 Institute of Information Transmission Problems of the Russian Academy of Sciences, Moscow, Russia, 

4 Moscow Institute of Physics and Technology, Dolgoprudny If 1700, Russia 
5 Universite Paris-Sud/CNRS, LPTMS, UMR8626, Bat. 100, 91405 Orsay, France, 

6 P.N. Lebedev Physical Institute of the Russian Academy of Sciences, 119991, Moscow, Russia 

We have considered an equilibrium ensemble of large Erdos-Renyi topological random networks 
with fixed vertex degree and two types of vertices, black and white, prepared randomly with the 
bond connection probability, p. The network energy is a sum of all unicolor triples (either black 
or white), weighted with chemical potential of triples, p. Minimizing the system energy, we see 
for some positive p formation of two predominantly unicolor clusters, linked by a ’’string” of Nbw 
black-white bonds. We have demonstrated that the system exhibits critical behavior manifested in 
emergence of a wide plateau on the Nbw (/i)-curve, which is relevant to a spinodal decomposition in 
1st order phase transitions. In terms of a string theory, the plateau formation can be interpreted as 
an entanglement between baby-universes in 2D gravity. We have conjectured that observed classical 
phenomenon can be considered as a toy model for the chiral condensate formation in quantum 
chromodynamics. 


I. INTRODUCTION 


During last two decades it has been recognized that 
Landau classification of phase transitions does not cover 
all patterns of symmetry breaking. The fractional quan¬ 
tum Hall effect and the Levin-Wen string networks pQ 
provide well known examples of topological order for so- 
called ’’gapped states”. Corresponding quantum phase 
transitions exist at T = 0 and do not have any local 
order parameters. Some symmetry aspects of 2D topo¬ 
logical phases with Z 2 symmetry have been discussed in 
[2J- Quantities describing topological phases are: i) the 
degeneracy of the ground state (which depends on the 
topology of the system), and ii) the holonomy of the 
non-Abelian Berry connection. Another way to identify 
’’topologically ordered” states is to use the topological 
entanglement entropy [3]. In any cases, the topological 
order is manifested in emergence of long-range correla¬ 
tions in the system. The interplay between some ’’local 
rules” and ” global constraints”, influence the behavior of 
the whole system and is responsible for the long-range 
order. 

Colored random networks have become an ubiquitous 
new paradigm for wide range of physical and social phe¬ 
nomena in distributed systems, spread from producers- 
consumers relations to string theory, and from budding 
in lipid membranes to baby-universes creation in cosmol¬ 
ogy. Here we study random black-white vertex networks 
with Levin-Wen type of the Hamiltonian [I], evolving by 
bond reconnections and tending to increase the number 
of unicolor triples under vertex degree conservation. The 
’’local rule” is the conservation of the vertex degree in 
all nodes of the network, while the ’’global constraint” 
is the attempt to maximize the number of triples of the 
same color. As the energy of triples increases, the net¬ 
work breaks into two predominantly unicolor (black and 


white) clusters, connected by a bunch of links between 
black and white vertices. Unexpectedly, this bunch is 
stable in a wide interval of triples energy. We conjecture 
that this novel statistical phenomenon has origin relevant 
to quantum entanglement. 

Ensembles of random Erdos-Renyi topological graphs 
(networks) provide an efficient laboratory for testing vari¬ 
ous collective phenomena in statistical physics of complex 
systems, being also tightly linked to random matrix the¬ 
ory. Besides investigating ’’typical” statistical properties 
of networks, like vertex degree distribution, clustering 
coefficients, ’’small world” structure and spectra of ad¬ 
jacency matrices, last two decades have been marked by 
rapidly growing interest in more refined graph character¬ 
istics, such, as for example, distribution of triadic motifs 
in oriented networks (small subgraphs involving different 
triads of vertices). 

Triadic interactions, being the simplest interactions be¬ 
yond the free-held theory, play crucial role in network 
statistics. Apparently, just presence of triadic interac¬ 
tions is responsible for emergence of phase transitions in 
complex distributed systems. First signature of a phase 
transition in a random network, known as Strauss clus¬ 
tering model 0, has been treated by RMT in [5] and 
identified with the 1st order phase transition in frame¬ 
works of mean-held cavity-like approach in jBj. Another 
example of phase transition is connected with the triadic 
motifs pattern formation, known as ’’motifs superfami¬ 
lies” , in real evolutionary networks J7]. The last prob¬ 
lem has been theoretically analyzed in [3j, where it was 
conjectured that stable motif prohles constituting super¬ 
families [71 may correspond to stability islands associated 
with localized states in the space of motifs. The local¬ 
ization of states occurs as a 1st order transition if the 
chemical potential associated with triadic motifs energy 
is large enough. The critical behavior of triadic motifs 
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concentration, as a function of the chemical potential, 
has been studied in details in J8]. 

Specifically, in this work we investigate the depen¬ 
dence of an equilibrium number of links, Nb w , connecting 
’’black” and ’’white” vertices in networks with fixed ver¬ 
tex degree, p 7 as the function of the chemical potential, 
p 7 controlling the number of unicolor vertex triples (all 
three vertices are black, or all three are white). The at¬ 
tempt to minimize the energy of unicolor triples leads to 
almost immediate color separation, when at any infinites¬ 
imal positive value of p two clusters with opposite colors 
(’’^ 2 -charges”) are formed. More striking phenomenon 
is observed for large enough / jl: the Nb w (p)-curve has a 
wide plateau separating two Arrenius-type dependencies 
at low and high p. Such a plateau is never seen (for 
any p) if the energy of the system is the sum of unicolor 
vertex pairs (but not triples), or if the vertex degree is 
not conserved. We suggest two parallel interpretation of 
the phenomenon observed in our simulations: ’’classical” 
(in terms of classical 1st order transitions), and ’’quan¬ 
tum” (in terms of spins interacting with 2D topological 
gravity). Describing the system in terms of classical sta¬ 
tistical physics, we interpret the plateau emergence as the 
spinodal decomposition similar to the vapor-liquid first 
order phase transition. 

On the other hand, the ’’quantum gravity language” 
implies the realization of the network as the triangula¬ 
tion of 2D fluctuating surface. Since the network does 
not involve metric, one deals with purely topological ver¬ 
sion of 2D gravity where the chemical potential of uni¬ 
color triples, p, can be understood as a 2D cosmological 
constant. The spin (i.e. color) separation phenomenon 
is interpreted as an anomalous transport on the string 
worldsheet due to the axial anomaly in the gravitational 
field. In terms of 2D gravity, it is conjectured that the 
plateau formation occurs as the specific ’’entanglement” 
of two baby universes. We propose the synchronization 
mechanism of phase transitions in each universe and sug¬ 
gest that their entanglement results in a sort of a phase 
coexistence phenomenon. To describe the phase coexis¬ 
tence properly, one could try to find gravity counterparts 
of thermodynamic variables. This is a subtle issue and 
we rely on recent interpretation of (P, V) variables in 
cosmological terms |9j deling with Van-der-Waals ther¬ 
modynamics of charged or rotating black holes m- 

We also mention useful analogy with the chiral sym¬ 
metry breaking in the instanton-antiinstanton ensemble 
in quantum chronrodynamics (QCD) described by ran¬ 
dom matrix model. In this case case, the system at small 
/i and low temperature is in the phase with the spon¬ 
taneously broken chiral symmetry characterized by the 
nonvanishing chiral condensate, ^’FlI/). Staying at T = 0 
and increasing the chemical potential, p, we force the 
system to undergo a quantum phase transition at some 
p c . The density of chiral condensate is proportional to 
the number of instanton-antiinstanton connections, and 
the chiral symmetry restoration signifies emergence of 


another condensate, ('I'd/). We conjecture that our two- 
color network could be a toy model for the formation of 
such condensate. 


II. THE MODEL AND KEY OBSERVATIONS 

The setup of the model is as follows. Take N vertices 
and label them by integers 1,..., N. Vertices i = 1,..., cN 
(0 < c < 1) are ’’black”, i.e. are associated with Ising 
spins <Ji = —1, while vertices i = cN+1, ...,N are ’’white” 
and carry cq = +1. The initial realization of Erdos- 
Renyi network is prepared by connecting any randomly 
taken pair of vertices with the probability p , {p > A, 
i.e. p is above the percolation threshold). Then, one ran¬ 
domly chooses two arbitrary links, say, between vertices 
A and B ( A-B) and between C and D ( C-D ), and re¬ 
connect them, getting new links A-C and B-D. Such 
reconnection conserves the vertex degree HU• Now one 
applies the standard Metropolis algorithm with the fol¬ 
lowing rules: i) if after the reconnection the number of 
unicolor vertex triples is increased, a move is accepted, 
ii) if the number of unicolor vertex triples is decreased by 
AW r ipi, or remains unchanged, a move is accepted with 
the probability g-MAAftripi _ w e have run the Metropo¬ 
lis algorithm repeatedly for large set of randomly chosen 
pairs of links, until it converges. For one-color networks it 
was proven [12] that such Metropolis algorithm converges 
to the Gibbs measure e AtJVtripl in the equilibrium ensemble 
of random undirected Erdos-Renyi networks with fixed 
vertex degree. For two-color networks the convergence 
is not proven rigorously, however extensive numeric tests 
do not show any pathologies. 

The obtained results are as follows. The dependence 
of number of links connecting black and white vertices, 
Nbw{p)i demonstrates for small and large p the Arrenius 
activation kinetics, N bw (p) ~ where /3(p) is some 

connectivity-dependent constant. However, the function 
Nb w (p) develops a plateau in a wide region of p, mani¬ 
festing a kind of a critical behavior. The corresponding 
plot for the network of N = 500 vertices is shown in the 
Figjl] for p = 0.15. Analyzing the clustering structure 
of the network for different p, one sees that at small p 
(p ~ ^), the system splits for any p into two clus¬ 
ters of predominantly black and white vertices, meaning 
the spontaneous Z 2 -symmetry breaking (the clusters are 
allocated by preferential number of links entering ’’in¬ 
side” it, rather than going ’’outside”). So, the plateau is 
developed in the Z 2 -broken phase. These results are re¬ 
producible for various values of c (where c is the fraction 
of black vertices). The same behavior is seen for random 
regular networks, which have fixed vertex degree, d, in 
all network nodes. 

Remind that we consider very dense networks, p = 
0(1), far above the percolation threshold, p = Cer¬ 
tainly, even in this regime the entire network contains 
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some number of disjoint clusters, exponentially sup¬ 
pressed in their sizes. In colored network, weights of 
different rewirings are color-dependent and could affect 
the distribution in cluster sizes. We believe that for dense 
networks (considered here) effect of color-dependent re¬ 
distribution in cluster sizes is negligible against the back¬ 
ground of exponentially suppressed overall number of dis¬ 
joint components. However, near the percolation thresh¬ 
old this effect might be visible and definitely deserves 
special attention. 

|i=0.0 |r=0.025 (.1=0.2 |r=0.3 (i=0.35 (i=0.4 



Figure 1: Dependence of equilibrium number Nbw of black- 
white bonds on the chemical potential, p, per each unicolor 
triple for N = 500-vertex random network of the probability 
p = 0.15 with equal composition of black and white nodes; 
the typical network topologies for some values of p along a 
Nbw(p)-c urve are shown above the Nbw(p)- curve. 

The following two key result are the subjects of our 
discussion and interpretation in terms of lattice gases and 
string theory: 

• Immediate spontaneous formation of black and 

white clusters (breaking Z 2 -symmetry) at p ~ ; 

• Emergence of a wide plateau in Wu) (/^-dependence 


III. SYMMETRY BREAKING 


system the interaction energy of one 6-node with m ran¬ 
domly chosen other nodes is Ub = jj(Nbfob + N w fo w ), 
where fob and fo w are the b-b and b-w interaction ener¬ 
gies. Respectively, Ub = ^(N b fob + N 2 <t> w b)- The free 
energy of the system reads F = \(NbUb + N w u w ) — 
kT (N b In ^ + N w In Af ). Defining </> = fo b + fo w and 

A = fob — fow , passing to concentrations, Cb lW = 

(Cb + c w = 1), and introducing the asymmetry of the net¬ 
work composition, r/ = Cb — c w , we can write the network 
free energy, / = (i, as follows: 


f(v) = fo+ ^ [(! - V) ln(l - V) 


+ (1 + rj) ln(l + r})] (1) 


l (f> + kT In 2. 


where /o = 4 ■ 

The system described by the thermodynamic potential 
([I]), experiences the standard phase transition with the 
stationary states, r] c , determined by the equations = 

0 and > 0. Supposing that \rj\ <C 1, we can 

expand f(rj) up to the 4th term and get the equation 


( toA \ kT o 

{— ~ hT )i~ Y’’ =° 


( 2 ) 


having a single solution rj c = 0 at m < ^ and exhibiting 
spontaneous symmetry breaking into two phases r/^ ,b = 

±v/ 3 (tY - !) atm> ^T. 

The mean-field arguments for accounting a plateau, 
are provided in the next section in terms of lattice gas 
statistics, random matrix and string theories. Being ap¬ 
plied to networks, the random matrix approach means 
the discretization of a two-dimensional string worldslieet 
embedded in some D-dimensional target space. Since 
the colored (black and white) vertices can be described 
by Ising spin variables, our system is some specific model 
of a matter (Ising spins) coupled to the gravity (fluctuat¬ 
ing surface over ensemble of allowed reconnections in the 
random network). In the matrix model framework, which 
parallels the mean-field analysis, such kind of systems has 
been considered in El. however with some important dif¬ 
ference: in our system Ising spins are quenched, while in 
matrix approaches they are typically annealed. 


A. Statistical interpretation 

First of all, it is worth pointing out that the Z 2 - 
symmetry breaking and the plateau formation are dif¬ 
ferent physical effects. Separation of clusters is typi¬ 
cal phenomenon which can be understood using naive 
mean-field arguments. Let Nb and N w be the number 
of black and white nodes (Nb + N w = N). In a mixed 


B. Stringy interpretation 

Nowadays the string theory plays two different roles 
in the fundamental physics. First of all, it is a model 
of quantum gravity since the spectrum of closed strings 
involves massless graviton. Secondly, the string can be 
regarded as the ’’probe”, and the physics on a 2D string 
worldsheet should reflect all phenomena happen in entire 
target space, in which this worldsheet is embedded. This 











4 


is a key idea behind the gauge/string duality. Having 
in mind these roles of string theory and considering a 
network as a discretized string worldsheet, we rephrase 
in stringy terms the phenomenon observed in the black- 
white network. 

The spontaneous symmetry breaking and the black- 
white cluster separation can be interpreted as occurrence 
of the anomalous transport on the string worldsheet. The 
anomalous transport emerges due to the quantum non¬ 
conservation of classically conserved currents in the ex¬ 
ternal gauge or gravitational fields. The spin separation 
phenomenon in our problem means that non-vanishing 
axial current takes place at the worldsheet when the spins 
in sub-ensembles S = {sigma = +1} and S = {a = — 1} 
move in opposite directions. Since spins are quenched 
and do not fluctuate, the current is induced purely by 
fluctuations of the 2D surface. The network is a topo¬ 
logical object, hence it is natural to describe it as the 
topological 2D gravity of Jackiw-Teitelboim type which 
involves the cosmological constant A 

L = J d 2 x^/g$(R + A), (3) 

where $ is the scalar dilaton held and R is the Ricci cur¬ 
vature of the two-dimensional metric g. This Lagrangian 
can be written as the topological 2D Yang-Mills theory 
with SL(2,R) gauge group (see [T3] for review). 

It is known that in the pure 2D gravity the total num¬ 
ber of triangles measures the total area, hence the cor¬ 
responding chemical potential, g, of triples can be in¬ 
terpreted as an effective 2D cosmological constant. At 
the equation of motion we have R = —A. There is an 
axial current due to the chiral anomaly in an external 
gravitational held 

d„4 = R ( 4 ) 

which for the equation of motion equals to the cosmo¬ 
logical constant, /i. Thus, for the axial current itself in 
terms of the spin-connection lo v we have 

4 = oc /a (5) 

Eq. ([5]) demonstrates the emergence of spin current for 
any nonzero chemical potential of triads, /x, which is the 
microscopic mechanism of black and white nodes separa¬ 
tion. 

The Z 2 spontaneous symmetry breaking in colored dy¬ 
namic network we describe using mean-held statistical 
arguments, as well as from more involved viewpoint, in 
terms of fluctuating surfaces. In what follows we pass 
along the edge of these two different approaches. 

It should be pointed out that spin separation phe¬ 
nomenon is known in other systems as well. At low tem¬ 
perature the 2D Ising model exhibits spontaneous magne¬ 
tization, developing the coherent spin domains separated 
by the domain walls. The Z 2 symmetry is spontaneously 


broken in this phase. In the case of fluctuating Ising 
spins coupled to 2D gravity, the spontaneous creation of 
’’baby universes” takes place, and it is known that the 
spin direction in the daughter universe is always oppo¬ 
site to the spin directions in the parent one (see [I^ for 
review). In our case we see the formation of clusters 
with opposite spin signs due to dynamic rearrangement 
of bonds only (since the spins are quenched). Similar 
spin separation phenomenon occurs in 2D matter with 
non-vanishing chemical potential (161 and in 4D matter 
in an external magnetic field (see m for review). The 
chirality separation emerges in some lattice QCD systems 
when the sheets of opposite chiralities become connected 
by stringy ’’skeletons” [18] . 


IV. PLATEAU FORMATION 

The most intriguing question concerns the interpreta¬ 
tion of the plateau formation, since it seems to be a phe¬ 
nomenon with the signature of a quantum phase transi¬ 
tion. In contrast to the symmetry breaking, the plateau 
formation is an essentially collective effect which disap¬ 
pears for networks with nonconserved vertex degree. 

The relevant image, which highlights the basic features 
of this phenomenon, is the famous puzzle ’’the game of 
fifteen”. The analogy goes as follows. When the concen¬ 
tration of graph links is small, it is always possible to 
minimize the system energy (the sum of unicolor triples) 
almost locally. This resembles the initial stages of order¬ 
ing in the puzzle, when one takes care of single-dice place¬ 
ments only. However, our network gets highly frustrated 
due to the vertex conservation condition: in order to find 
the way to place some new ” good” link which minimizes 
the unicolor triple energy, it might be necessary to pass 
over high potential barrier and remove other ’’good” links 
which have been already placed before. Thus, the en¬ 
ergy minimization involves unfavorable sequences of link 
permutations until all necessary constraints are satisfied. 
This is exactly what happens on late stages of ordering 
in the ” game of fifteen”. 


A. Statistical interpretation 

To proceed with combinatorial interpretation, repre¬ 
sent the graph by the symmetric adjacency matrix Rj 
and split it in four sectors with corresponding numbers 
of particles (links) N bb , N bw , N wb , N ww : 


N bb in [••] for [1 < {?’, j} < Nc] _ 

N ww in [00] for [Nc+ 1 < {i, j} < N} _ 

N bw in [*o] for [Nc+l<i<N, 1 < j < Nc] 
N wb in [o*] for [1 < i < Nc , Nc + 1 < j < V] 
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The matrix is symmetric, namely Nbw = N w b- The 
unicolor triple is any pair of particles in one row or in one 
column in the sectors [••] or [oo]. The vertex degree con¬ 
servation implies the condition = ^j°\ where 

Jj 0 ' 1 (J = 1, ■■■! N) are the initial vertex degrees, fixed by 
the network preparation conditions. Hence the simplest 
rearrangements of graph links should involve pairs of par¬ 
ticles (links) shown in the Fig|2] As we shall see, this 
figure helps to understand the appearance of the plateau 
on the Nb w (^)-curve depicted in the Fig[l] 
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Figure 2: (a) Pulling particles (bonds) 1 and 2 (and their 
symmetric counterparts 1' and 2') from the sectors [«o] and 
[o»] into the sectors [••] and [oo]. The rearrangement results 
in two newborn triples of graph vertices (pairs of particles); 
(b) Bond rearrangement which conserves number of particles 
in all sectors, however increases the number of triples in the 
sector [••]. 


The physics behind the Nbw (/^-dependence is as fol¬ 
lows. When the chemical potential, /t, of unicolor triples 
is increasing, the formation of aligned pairs of particles 
(in rows or in columns) in the sectors [••] and [oo] be¬ 
comes favorable. There are two mechanisms which in¬ 
crease the number of aligned pairs in these sectors, they 
are schematically shown in the Figs|2ji,b 


1. In the first (I) mechanism the exchange of particles 
(shown by arrows in the Fig{2^) 1^2 (and of their 
symmetric counterparts, V 2') pushes particles 
1 , V to [••] and particles 2,2' - to [oo] 

2. In the second mechanism (II) - see the Fig§ 3 
the exchange 1^2 (and V ^ 2') preserves the 
number of particles in all sectors. 


At small fj, (at the beginning of the Nb w (p)- curve) 
both mechanisms (I and II) of unicolor triples energy 
increasing are available. However, as the number of 
particles in the sectors [••] and [oo] exceeds some crit¬ 
ical value, the entropic loss due accumulation of parti¬ 
cle in these sectors forces the 1st order vapor-liquid-like 
phase transition in the sectors [••] and [oo]. Note that 
since the initial concentrations, /(N 2 c 2 ) and 

Cwh = Nwl/(N 2 ( 1 — c) 2 ), of particles in sectors [••] and 
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Figure 3: Number of unicolor triples, AAripi as a function of 
/i in the ensemble of networks with N = 500 and p = 0.15. 


[oo] do not coincide in general (they are fixed by a ran¬ 
dom network preparation), the difference 


r (°) _ r (o) 

L 'bb ^ww 


(Np) 


- 1/2 


( 6 ) 


desynchronizes the vapor-liquid transitions in [••] and 
[oo]. So, the beginning of the plateau on the N bw (p)~ 
curve we connect with the first occurrence of independent 
vapor-liquid transitions in [••] or in [oo]. The plateau 
itself we identify with the synchronization of these two 
vapor-liquid transitions in [••] and in [oo]. At the plateau 
the mechanism I gets suppressed, because capturing of 
particles by sectors [••] and [oo] becomes entropically 
unfavorable. 


The plateau on the N bw (/r)-curve occurs simultane¬ 
ously, as shown in the Fig[3j with the plateau on the 
curve -/V tr ipi (/t), where 7Vt r i P i is the number of unicolor 
triples. The plateau on the A^dpi (/i)-curve means the 
topological quench of each unicolor subnetwork in [••] 
and in [oo]. This behavior is a signature of a spinodal 
decomposition typical for the vapor-liquid phase transi¬ 
tions. The mechanism which competes with increase of 
energy of triples in the plateau region is just the mecha¬ 
nism II of particle exchange: particle pairings, as shown 
in the Figj2j}, do not change the numbers of particles in 
all sectors, however act against the entropy in the sectors 
[•o] and [o*] since trapping of particles like 2 and 2! is 
entropically unfavorable. 

Collecting all together, we are led to the following 
conjecture. The plateau begins at the value /if , eg inde¬ 
pendently in subnetworks [••] and [oo] as the individual 
vapor-gas phase transitions. At the plateau the interac¬ 
tion between the sectors [••] and [oo] is ensured by the 
process II, particle trapping in [••] and [oo] works against 
the entropy in [*o] and [o*], and we deal with two effec¬ 
tively coupled subsystems [••] and [oo], both exhibiting 
Is order transition and living in a joined metastable state. 
The plateau is finished at the value p en d, at which both 
sectors [••] and [oo] fall down in liquid states, and de¬ 
couple. At /t > Pend the entropic contributions in all 
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sectors become negligible and the Arrenius-type activa¬ 
tion kinetics gets restored. 

Note, that just above the percolation threshold, i.e. for 
p > jj the standard Erdos-Renyi network, besides the 
giant component, contains many disconnected clusters 
with the size distribution c(k) ~ fc -5 / 2 , where c(k) is the 
concentration of clusters of k vertices. For black-white 
networks the distribution c(k) should be modified since 
weights of particular clusters become color-dependent. 
For all p clusters contain only nodes of the same color, 
and all black-white links are external, connecting oppo¬ 
sitely colored clusters. The very existence of the plateau 
is insensitive to the number of disjoint unicolored clusters 
in the network, however the plateau disappears for quite 
small p. 


B. 2D gravity/stringy interpretation 

In terms of 2D gravity, the emergence of the plateau 
can be linked with a ’’baby universe” formation. It has 
been proven in m ns that critical phenomena in 2D 
gravity deal with the formation of a newborn (’’daugh¬ 
ter” ) bubble connected by the ” neck” to the parent uni¬ 
verse. The free energy of the system consists of two com¬ 
peting contributions: the mean value of the total ’’neck 
diameter”, and the spontaneous curvature of the new¬ 
born (baby) universe of a given size. The same mecha¬ 
nism lies behind the ’’budding of vesicles” in liquid mem¬ 
branes [20]. In terms of the black-white network, the 
color separation can be interpreted as the baby universe 
formation, whose size depends on the ratio /Nww, 

where N^' (or Nww) are the initial numbers of black- 
black (or white-white) nodes, while the interaction of 
triples plays the role of a curvature contribution to the 
free energy (Exed by the cosmological constant, p). 

The plateau formation, in terms of 2D gravity, signals 
that the neck size (or sum of all neck sizes for multi¬ 
cluster network) connecting the parent and baby uni¬ 
verses, does not depend in some interval on the bulk 
cosmological constant, p. The pure 2D gravity under¬ 
goes the phase transition at some value p c (see j21j for 
review), therefore upon the color separation, we have two 
copies of 2D gravities connected by the neck. Each grav¬ 
ity could be in one of two phases. Since we consider 
the quantum geometry, there is no reason for initial syn¬ 
chronization of phase transitions in these two universes. 
The value of the cosmological constant, at which the first 
transition happens, is presumably the point at which the 
plateau gets started. Two universes are entangled, which 
in terms of each 2D gravity can be interpreted as the in¬ 
sertion of the macroscopic loop operator. 

So, the phase transition in one (say, ’’black”) universe 
happens spontaneously, however due to the quantum en¬ 
tanglement of two universes ([••] and [oo] clusters), the 


transition in the second (’’white”) universe is induced 
by the vacuum expectation value (vev) of the macro¬ 
scopic loop operator. The plateau corresponds to the 
coexistence of two phases in the 2D gravity and when 
the second universe undergoes the induced phase tran¬ 
sition, they get separated and the plateau terminates. 
Note that we consider the Euclidean version of the 2D 
gravity, when the baby universes can be produced. 

If we would deal with fluctuating Ising spins coupled 
to 2D gravity, the results could be borrowed from the 
literature. The case when spins live on faces, has been 
solved in na, while the dual Ising model for spins living 
on vertices, was treated in [22]. In the last case the model 
has been reduced to the 0(1) loop model. In our system 
we have spins on the vertices, hence the dual Ising model 
1221 is more relevant. It admits the two-matrix represen¬ 
tation [22] 

Z = J d XdYe Tlx2+ ^ x3+TrY2+ ^ XY2 (7) 

where X is the adjacency matrix both for black-black and 
white-white links, while Y is the matrix of black-white 
connections. Besides, we have two additional ingredients 
compared to [22] : the spins are quenched, and the graph 
vertex degree is conserved. The degree conservation can 
introduced into the matrix model by the linear term with 
the Lagrangian multiplier TrMX, playing a role of an 
external magnetic field similar to [23] . At Y = 0 we are 
back to the pure matrix model describing the 2D gravity 
with the cosmological constant p. 

The plateau formation has an analogy with the chiral 
symmetry breaking in QCD with non-vanishing baryonic 
chemical potential, p. Remind few key QCD properties 
relevant to our system. At T = 0 the QCD is in the 
confinement phase with chiral condensate, ('FT) ^ 0. 
The popular model of the QCD vacuum is the instanton- 
antiinstanton (II) liquid, represented by an ensemble of 
interacting point-like objects of two types. The value 
of chiral condensate is presumably expressed in terms of 
number of connected instanton-antiinstanton pairs [24] . 
which is the number of black-white links in the network 
language. The number of fermionic zero modes at the 
instanton is an unchanged topological number, mean¬ 
ing fixed vertex degree in network terms. In the II- 
ensemble zero modes at individual instantons get col¬ 
lectivized. The phase diagram in QCD is quite rich in 
the ’’temperature-chemical potential” plane. At small p 
the system is in the confinement phase with broken chi¬ 
ral symmetry. By increasing p : one reaches mixed phase 
with two coexisting condensates, (TT) and (TT). Such 
a mixed phase exists in some region in p , in which the 
string tension does not depend on the density, being the 
analogue of the network plateau. Further growth of p 
forces the system to fall down into the so-called color 
superconductor phase, with chiral symmetry restoring 
and QCD string disappearing. In the QCD color super¬ 
conductor phase, which is believed to occur in neutron 
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stars, condensation of quark-quark state happens similar 
to the Cooper pairing in conventional superconductors. 
The phenomenon of chiral symmetry restoration at large 
/r apparently can be described by the effective random 
matrix model (see [55] for review). We conjecture that 
the two-color network picture could be regarded as a toy 
model for zero-mode collectivization and can be related 
with holographic representation of the chiral symmetry 
breaking proposed in 5(Jj. 


V. DISCUSSION 

In this paper we have described a novel critical phe¬ 
nomenon in statistics of two-color random network - the 
emergence of a wide plateau in the concentration of bonds 
linking black and white nodes. The key condition for the 
plateau development is the graph vertex degree conser¬ 
vation, the very natural condition for any topological en¬ 
semble. Physics behind the phenomenon can be thought 
of as the particular mechanism of bonds collectivization 
in Levin-Wen-type topological Hamiltonians. It has evi¬ 
dent parallels with the effect of fermionic zero modes col¬ 
lectivization, emerging from localized solutions of Dirac 
equation in the instanton-antiinstanton QCD ensemble. 
Remind, that this effect results in chiral symmetry break¬ 
ing at small condensate density and its restoration at 
large density. It would be highly desirable to find an 
appropriate matrix model description of phenomenon, 
similar to the QCD chiral matrix model describing the 
spectrum of the Dirac operator in the II ensemble. 

One could consider ensembles of topological defects of 
different nature like instantons, monopoles or vertices. 
The effect of the plateau formation due to the zero modes 
collectivization is expected to be universal for them. Re¬ 
cently, formation of the condensate in the ’’colorless” 
instanton ensemble without antiinstantons was investi¬ 
gated where it has been found that microscopic de¬ 
scription of condensate involves refined knot invariants. 
The critical behavior in such ’’colorless” topological en¬ 
sembles inspires the conjecture that the knot recogni¬ 
tion by the topological invariant is different in differ¬ 
ent phases. It would be very interesting to understand 
whether the relation with the knot invariants can be 
found for colored topological ensembles. 


In this work we follow the probe approach for the 
string, however one could question whether the critical 
plateau formation would have stringy meaning if consid¬ 
ering the colored network as a quantum gravity model. 
The idea to use colored network in this way has been 
recently suggested in [28]. However our model differs 
from the one proposed in 1281 in two crucial respects: i) 
we consider conserved vertex degree, and ii) in our case, 
contrary to (28j . the flow increases the number of triads, 
since the signs of cosmological constants in our model 
and in [2S] are different. 

The novel criticality found in colored networks seems 
to be quite general phenomenon and could have the prac¬ 
tical applications in real life networks where there are 
sets of separated communities (black and white nodes). 
The criticality takes place even in highly asymmetric 
black-white networks (i.e. for the parameter c lying in 
the region about [0.1,0.9]). As an immediate example 
we could mention social networks. In this framework, 
the phenomenon of the plateau formation sounds as fol¬ 
lows. Consider two social communities, and assume that 
each community desires to increase the number of tri¬ 
adic ’’own connections” (simplest ’’own cliques”) mea¬ 
sured with some weight (chemical potential of triads). 
We predict in such networks existence of the ’’stability 
plateau” (i.e. presence of unavoidable communication 
between communities), at which the number of links be¬ 
tween different communities is insensitive to the weight 
of ’’own connections”. 
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